PROCEEDINGS 

0F SCIENCE 



Electromagnetic corrections to light hadron masses 



h A. Portelli*', S. Durr ad , Z. Fodor acrf , J. Frison fe , C. Holbling S.D. Katz QC , 

S. Kr\eg ad , T. Kurth a , L. Lellouch 6 , T. Lippert K.K. Szabo" and A. Ramos z 
CN (Budapest-Marseille-Wuppertal Collaboration) 

"Bergische Universitdt Wuppertal, Gauflstr. 20, D-42097 Wuppertal, Germany 
b Centre de Physique Theorique, Luminy, Case 907, F-13288 Marseille cedex 9, France t 
c Institute for Theoretical Physics, Eotvos University, H-l 117 Budapest, Hungary 
d Jiilich Supercomputing Centre, Forschungszentrum Jiilich, D-52425 Jiilich, Germany 



(N 



c3 



> 



X 



E-mail : portelliScpt . univ-mrs . f r 



Oh 
D 

At the precision reached in current lattice QCD calculations, electromagnetic effects are becoming 
numerically relevant. We will present preliminary results for electromagnetic corrections to light 
hadron masses, based on simulations in which a U(l) degree of freedom is superimposed on 



ON Nf = 2 + 1 QCD configurations from the BMW collaboration. 
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1. Motivation 

Isospin is a near symmetry of the hadron spectrum because the strong interaction only distin- 
guishes quark flavors through their masses and the mass difference between up and down quarks is 
small. This symmetry is broken by : 

(i) the mass difference m u — m^ (mass isospin breaking) 

(ii) the difference in the charge of the u and the d quark (electromagnetic isospin breaking) 
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5 05+ - 75 



These effects are expected to be at the percent level. The size of mass breaking is the mass differ- 
ence m u — m c i relatively to a typical QCD scale Aqcd and the order of electromagnetic breaking is 
the fine structure constant at zero momentum ce = ^ ~ 

These effects imply mass splittings inside isospin multiplets. Although these effects are small, 
they have important consequences. For example, the fact that the neutron is heavier than the proton 
guarantees the stability of matter. Another interesting isospin breaking quantity is the absolute 
correction to Dashen 's theorem : 



A A D = AemMk - A EM M^ 



(1.1) 



where 



A EM M| = (Mp+ -Mp„)m„= 



-111,1 



(1.2) 



is the electromagnetic squared mass splitting of the isospin multiplet P. One can also consider the 
dimensionless relative correction to Dashen 's theorem : 



AhmM|_ 1 
A EM M2 



(1.3) 



Dashen has shown in [2] than A&D = in the SU(3) chiral limit and that the leading corrections 
are 0(ccm s ,GC 2 ). The quantity A&D is interesting because it is very sensitive to the up and down 
quark masses. 

An overview of the results for the violations to Dashen's theorem collections is presented in 
Table 1 . These results are interesting for two reasons : on one hand, agreement between different 
calculations is poor and on the other all these numbers predict rather large corrections. We will 
investigate corrections to Dashen's theorem using a QCD+QED analysis on BMW collaboration 
QCD ensembles. A first step in this direction is to formulate electromagnetism on the lattice. 
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1230 
1300 ±400 

360 
1060 ±320 

1080 

1070 



526 
340 ±92 
1250 ±550 



ArD 



0.80 [3] 

1.02 ±0.30 [4] 

0.26 [5] 

0.84 ±0.24 [6] 

0.68 [7] 

0.85 [8] 

039 [9] 

0.30±0.08 [10] 

N/A [11] 



Table 1: Results for the violations to Dashen's theorem. Bold numbers are the results given by the authors, 
the others are deduced from information given in the corresponding paper. 

2. Electromagnetism on the lattice 

We impose periodic boundary condition in a finite volume, i.e. we will define Maxwell's theory 
on a 4-torus T 4 . In this context, the Maxwell-Gauss equations : 



£ duFuv = jv with F[i V = <9 M A V - d v A^ 
imposes global electric neutrality : 

Gtotai = / d 3 xj'o(x) = T I d 3 xd k F k0 (x) = 

JT 3 £fj JT 3 



(2.1) 



(2.2) 



One way to circumvent this global neutrality is to add a finite-volume, constant background current 
density as a multiplier to the Lagrangian of the theory : 



(2.3) 



M=0 



where is the length of direction /I in T 4 and V = L0Z4L2L3 is the 4- volume of T 4 . The Maxwell- 
Gauss equations become : 

1 

V~ 



y, d^F^y — j v — —L v c v 



(2.4) 



Minimizing the action with respect to the Lagrange multiplier c yields a new Euler-Lagrange equa- 
tion : 

44„„ f^_ n (25) 



A M (0)= / d 4 xA M (*)=0 

JT 4 



where the hat symbol denotes Fourier transform. 

We formulate now electromagnetism on the lattice in its non-compact form (following [9]). 
This is done by using forward finite differences in the Lagrangian (2.3) : 



d„/to = /(*+£) -/to 



(2.6) 



3 



Electromagnetic corrections to light hadron masses 



A. Portelli 



This formulation has a major advantage : the photon does not self-interact. The inconvenience is 
that gauge fixing becomes compulsory. We chose the following gauge conditions : 



Vpo^O, A (po,0) 
Vp, Ll=\PkAk(p) 








with pn 



2 sin — 
V 2 



(2.7) 



This choice has the practical advantage that it allows to solve explicitely the gauge constraints 
during field generation (cf. appendix of [10]). 

The action S associated with the above Lagrangian is a quadratic form, so generating random 
fields according to the law DA exp(— S[A]) to compute the path integral of the theory is simply a 
multi-dimensional Gaussian random number drawing. This is done in momentum space, where 
the correlation matrix of the law is quite sparse, and then the position representation of the field 
is recovered using a Fast Fourier Transform algorithm. The associated U(l) compact field (links) 
ufi ED is then constructed in the following way : 



77QED 



(2.8) 



where q is an electric charge. 

To check the implementation of the field generation code and investigate volume dependence, 
we computed semi-analytically the expectation value of the plaquette operator P^ v associated with 
the field uj} ED in an infinite volume lattice to compare it with simulation values. The results are 
summarized in Figure 1. It appears that for lattice extents greater than 10 the plaquette does not 
suffer from finite- volume effects and hypercubic invariance is present even with our non symmetric 
gauge choice (2.7). 



-4r log((Pp„» with q = 1.0 on L 4 lattices 
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0.49 - 



0. 185 
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Figure 1: Logarithm of the U(l) plaquette vs. lattice extent L. The dashed blue line represents the semi- 
analytical, infinite-volume prediction. The green points correspond to plaquettes aligned along the time and 
a space directions and the red ones to plaquettes aligned along two perpendicular space directions. 
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3. Coupling to quarks 

Here we describe how to perform a QCD and quenched QED analysis by re-using previously 
generated SU(3) gauge configuration. We are using Nf = 2 + 1 QCD simulations with Liischer- 
Weisz gauge action, tree level 0(a)-improved Wilson fermions and two steps of HEX smearing 
{cf. [12, 13] for details). 

Ideally one would include the U(l) degrees of freedom directly in the Hybrid Monte-Carlo. 
But doing this to reproduce the wide range of parameters already explored in BMW collaboration 
ensembles is extremely expensive. If, in a first instance, one accepts to quench QED, then there is 
a simple way to re-use previously generated SU(3) fields. To compute the propagator of a quark q 
of charge Q q G {3,-5} proceed as follows : 

1. generate an electromagnetic potential as described in Section 2 

2. construct the associated U(l) field £/J ED = exp(iQ q eA^) 

3. using a previously generated SU(3) field U^ CD , construct the U(3) field 

4. compute the quark propagator inverting the Dirac -Wilson operator using as the gauge 
field in covariant derivatives 

Because Wilson fermions break chiral symmetry explicitly, they require an additive renormal- 
ization in the quark masses. This additive renormalization has already been accounted for in our 
pure QCD simulation. However, in the presence of an extra U ( 1 ) gauge interaction, we expect an 
additional additive renormalization term, that will have to be subtracted, of order 0(<2^") where 
a is the lattice spacing. While suppressed by a, this divergence may still be large because it is 
enhanced by a -1 . Moreover, this correction breaks the mass isospin symmetry of our Nf = 2 + 1 
simulations, since Q u / Qj. We have chosen to subtract this divergence by retuning the bare u and 
d quark masses in such a way that the renormalized masses are equal. Thus, we substract a quantity 
QqS from each bare quark mass and tune S to a value S c where isospin symmetry is recovered. This 
isospin tuning works very well, as shown in Figure 2. 

Everything is now ready to perform a QCD+QED analysis in the mass isospin limit to inves- 
tigate the violations to Dashen's theorem. 

4. Preliminary results 

Here are our first preliminary results, using the following subset of the BMW collaboration 
QCD ensembles : 
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Mass isospin tuning (32 x 16 3 , = 3.31, mJJ' = -0.085, mf"'° = -0.04) 




Figure 2: PCAC light quark masses vs. 8. The behaviour in 8 appears to be perfectly linear. One can sees 
the large contribution of the a -1 additive renormalization at 8 — 0, where no substraction is applied. 

This subset contains a single lattice spacing estimated to be a ~ 0.116 fm (aT l ~ 1.7 GeV), one 
strange quark mass approximatively tuned to its physical value and four pion masses, from 420 MeV 
down to 200 MeV. Clover improvement is done only on SU(3) fields. 

Time correlators of n°, 7i + , K Q and K + are constructed by contracting quarks propagators 
computed as described in Section 3. The 71° correlator, is constructed by averaging the uu and 
d d pseudoscalar propagators, neglecting disconnected contributions. These contributions are sub- 
leading in a and m u — m c i. The mass of a pseudo-Goldstone meson is extracted fitting its time 
correlator to a hyperbolic cosine. The square of this mass is then extrapolated to the physical 7t + 
mass = 139.57018 MeV) using a Taylor expansion around (as done for [14]) : 



mXmI + )=m£' 



i+I>(m2 + -m; + 2 )* 



k=0 



(4.1) 



With this methodology, the previously described QCD ensembles and a linear extrapolation yield 
the preliminary results (quoted errors are statistical only) : 





= 134.5±l.lMeV 


Aem^ 


= 5.1 ±1.1 MeV 


AemMI 


= 1380 ±50 MeV 2 


M K + 


= 501.3 ±2.0 MeV 


M K o 


= 499.0 ±2.0 MeV 


&emM k 


= 2.2 ±0.2 MeV 




= 2200 ±180 MeV 2 


A A D 


= 830 ±180 MeV 2 


ArD 


= 0.60±0.14 
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5. Conclusion 

The preliminary results obtained in Section 4 validate our methodology and clearly encourage 
us to perform much more detailed calculations (lower quark masses, several lattice spacing, differ- 
ent volumes, etc. ) to gain control over systematic errors and to extend these calculations to other 
interesting hadronic observables. 
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